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Abstract 

We present a thermo-statistical model of popcorn production and 
propose a way to control the final size of the popcorn by monitoring 
only the chamber pressure. 
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How can one control the size of popcorn? The food industry may have 
ample reason to find ways to double or triple the size of popcorn, because 
large popcorn may be pleasing to the eyes, possibly energy saving in its mass 
production, and thus profitable to producers and consumers alike. The pur- 
pose of this note is to examine the physics of popcorn production based on 
a few thermodynamic principles and recommend ways to control the size of 
the final product. Before going into detail, let us first examine the underly- 
ing physics of popcorn production. First, there is moisture inside the corn, 
which, upon being heated in a chamber (e.g. a microwave oven), is subject 
to expansion. When the temperature of the chamber exceeds the boiling 
temperature, the moisture turns into a gas with a sudden expansion in its 
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volume. But because of the hard shell that surrounds the corn, the moisture 
is trapped inside and its pressure rises. When the gas pressure reaches the 
yield point of the shell, the shell is broken and the trapped gas undergoes a 
rapid adiabatic expansion, because the process is so fast that the gas does 
not have enough time to exchange heat and equilibrate itself with the envi- 
ronment. This adiabatic expansion stops when the gas pressure P reaches 
that of the chamber, P Q , which is, under normal circumstances, the atmo- 
spheric pressure. Since our primary interest here is to control the final size 
of the popcorn, and we know when the expansion stops, all we need do is 
write down the governing equations of the interface right after the rupture 
of the shell. Hence, this problem is reduced to that of interfacial instability 
and pattern formation, i.e. we are dealing with an interface advancing into 
the chamber. 

To write the equations, we recognize that the dynamic variable is obvi- 
ously the pressure P inside the corn at time t, which should obey the law of 
adiabatic expansion: 



where V(t) is the volume of the corn, and 7 is the ratio of the specific heats 
at constant pressure and at constant volume. For water vapor, 7 1.3 (1). 
The constant C Q is determined by the yield pressure, Py, of the shell and 
the initial volume of the corn, V Q) i.e. C Q = PyV~'. It is physical to assume 
that the advancing velocity of the interface, or the normal velocity, v n , is 
proportional to the pressure gradient: 



P(t)V(ty = C Q = const 



(1) 
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where k is a material constant. Eqs. (1) and (2) define the dynamics of the 
interface. 

We next examine the stability of the interface in a standard way. We first 
scale the time in such a way that k — 1. If the radius of the corn at time t is 
R(t), then its volume is simply given by V(t) = §7ri? 3 (t). The evolution of 
the sphere is then governed by Eq. (2): 

v n = dR(t)/dt = P-P = A/R 31 - P (3) 

where A = C G (^) 7 . Note that 7 is a non-integer, and thus Eq. (3) is highly 
nonlinear, and the exact solution is not available. However, the asymptotic 
form of the solution is easy to find. In the beginning, one may set P Q = 0, 
because P Q << P, and obtain an approximate solution, 

i?(t) 1+37 « A(l + 3 7 )t + 5 (4) 

where B = R l +Zl with R Q being the initial radius at t — 0. In the limit, 
t — > 00, the asymptotic solution is obtained by setting the left hand side of 
Eq. (3) to zero, i.e.: R(oo) = (A/ P ) 1 / 31 . So, R(t) initially increases as a 
power law, and then approaches R(oo). We have solved Eq. (3) numerically, 
and indeed checked this behavior. To examine stability of this solution, we 
imagine the sphere is perturbed slightly. Since any perturbations in three 
dimensions around the sphere can be expanded in spherical harmonics, Yi m , 
let 

r{t,e,4>)=R{t)+A lm {t)Y lm {9,4>) (5) 

If the amplitude A\ m grows in time, then the interface is unstable to the 
perturbation. Otherwise, it is stable. We should now obtain the corrections 
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to the pressure and volume caused by this perturbation. Define P(t) = 
P(t) + f(t) and V(t) —V(t) + g(t), where / and g will be functions of Ai m 
and Y lm . We first obtain g(t) in the case I ^ by direct integration of the 
volume: V(t) = § J> 3 (t)dfi = §7ri? 3 + RAf m . Hence, we find: g(t) = RA 2 lm 
and 

/(*) = (6) 

Now, in order to investigate the stability, we need a time dependent equation 
of motion for the amplitude A im . The normal velocity, v n , is given by: 

v n = dv/dt •n = P-P + f(t) (7) 

where r = (x, y, z) = r(t, 8, <fi) (sin 9 cos 0, sin sin 0, cos 0) with r(t, 9, <p) given 
by (6) and n is the normal vector to the surface, which is given by: 

n = (-dz/dx, -dz/dy, + (dz/dx) 2 + (dz/dy) 2 (8) 

Canceling R(t) = P — P D (Eq. (3)), and calculating the derivatives by chain 
rules yields: 



, m a v m ,, (g + A lm Y lm (6, <j>))((deY lm ) 2 sin 2 + (d^) 2 ) a2 

J it) = A lm Y lm (U,(p) 2i? 2 sin 2 A im + - 

(9) 

Since there is no angular dependence in f(t), we conclude from Eq. (9) that 
Aim = 0, namely all higher order harmonic perturbations of I ^ must be 
marginally stable. If we consider only the radial perturbation, then we may 
set, r(t) = R(t) + S(t). Then it is easy to show that the growth rate uo is 
negative: 

"^-^ + «<» do) 
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Hence, the interface is stable against the radial perturbation. In summary, 
we have shown at the level of the linear stability analysis, that perturbations 
of high order harmonics are marginally stable, and the lowest radial mode 
of I = decays exponentially. We now present a physical argument why 
the interface is indeed stable. Note that the volume correction due to the 
perturbation, g(t), is positive, while the correction to the pressure, f(t) is 
negative. Thus, when the perturbation arises, the inside pressure that drives 
the instability decreases, and the propagating speed decreases. Consequently, 
the instability is suppressed. Note that this problem is analogous to the 
problem of a solidifying interface advancing from a cold to a hot environment, 
which is always stable against perturbations. 

Since we have established that the interface is stable against infinitesimal 
perturbations, we are now in a position to determine the approximate size 
of the popcorn at a given chamber pressure P Q . Suppose the corn stops its 
expansion at a certain time tf, at which point P(tf) = P Q . We can easily 
find from Eq. (1) the maximum volume of the corn at tf\ 

V{t f ) = {C /P fh (11) 

Note that V(tf) is a function of the initial chamber pressure P Q , yield pressure 
Py and the initial volume, V , because C Q = PyVJ. For a given corn that 
has fixed V Q and Py, Eq. (11) enables us to control the size of the popcorn 
by monitoring only the chamber pressure P Q . Also, from Eq. (11), we may 
have a rough estimate of the yield pressure. Examining typical popcorn, we 
note that the radius increases at least by a factor 4, and thus Vf/V ~ 60. 
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Hence, 

Py/Po = (v f /v y « 200 

We note that the yield strength of the shell at the rupture point is about the 
same as or greater than polyethylene (LDPE) at room temperature (2). We 
now define the volume expansion rate T as the ratio of the final volume over 
the initial volume V Q : 

t(p ) = v(t f )/v = (c /p y/yv (12) 

If our goal is simply to increase or decrease the size of the popcorn relative 
to a given reference point, then the more relevant quantity is the ratio a = 
r(i^)/r(P ). From Eq. (12), we obtain: 

« = (Po/K) lh (13a) 

or equivalently, 

P' = a~"<P (13b) 

which is the central result of this paper. If we want to increase the size of 
the popcorn by a factor a, then we must reduce the pressure of the cham- 
ber by a factor a' . If a — 2, one must reduce the chamber pressure by a 
factor 2 13 2.5. However, we caution that the assumption of our model 
that the solid expansion of the corn via eversion closely follows the adiabatic 
expansion of the gas may be too simplistic, and thus the exponent, I/7, in 
(13a), which characterizes the expansion rate, may be somewhat smaller and 
presumably saturate beyond the critical pressure. We will report the experi- 
mental results in the near future. Another way of controlling the size of the 
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popcorn may be to increase the pressure P inside the corn and the cham- 
ber pressure P Q simultaneously, and let the chamber undergo a very rapid 
free expansion, when the temperature of the chamber reaches the popping 
temperature. During this free expansion, the chamber pressure decreases, 
while the pressure inside remains constant. When the difference, P — P Q , 
reaches the yield point, Py, the shell breaks, and the the corn undergoes the 
adiabatic expansion. By controlling the duration of the heating and the free 
expansion, one may control the size of the final product. This method may 
be in fact more promising than the previous one, because by not reducing 
the chamber pressure, it prevents moisture from being leaked to the chmaber 
during heating, which in turn increases the pressure of the corn, and hence 
the pressure gradient of the solid-gas interface during the expansion. 
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